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(54) Abstract Title 

A method, including modelling, for determining physical parameters of a sample 

(57) Method for faster determining certain physical parameters (capillary pressure curves and relative 
permeabilities for example) of a geologic sample initially saturated with a first fluid placed In a vessel and 
subjected to centrifugation in the presence of a second fluid. The equilibrium saturations are determined long 
before complete stabilization is reached, through analytic modelling involving exponential functions, best 
adjusting to an effective production curve, for each speed stage. A database obtained by means of various 
simulations is used to obtain the relative permeability curves corresponding to the centrifuged sample. The 
method can be applied for measurement of petrophysical characteristics. 

The method comprises measuring the amount of fluids displaced by centrifugation and establishing 
saturation curves, characterized in that it comprises, for each rotating speed, modelling a saturation curve best 
corresponding at each time interval to the curve portion already obtained. The figure shows a centrifugation 
device having a cup 13 for the porous rock sample, which is initially saturated with brine and placed in cup 13. 
The cup 13 is Initially filled with an electrically insulating fluid such as oil. 
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At least one drawing originally filed was informal and the print reproduced here is taken from a later filed formal copy. 
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OPTIMIZED METHOD FOR DETERMINING PHYSICAL PARAMETERS 
OF A SAMPLE SUBJECTED TO CENTRIFUGATION 

The present invention relates to a modelling method allowing to obtain more rapidly 
certain physical parameters of a geologic sample from e}q)erimental results obtained by 
subjecting it to centrifiigation. 

The method finds interesting applications in the field of petrophysics where 
knowledge of parameters such as the capillary pressure, the relative permeability curves 
of porous rocks, aUows reservoir engineers to characterize the flow conditions of 
hydrocarbons in a resarvoir, to make production predictions and to elaborate 
development strategies. 



A well-known method for rapidly determinii^ both the capillary pressure curves 
and the relative permeability curves of a porous rock sanq>le initially saturated with a 
first fluid essentially consists in subjecting it to centrifiigation in order to displace the 
fluid in place by means of another fiiud. 

15 Centrifiigation devices such as those described in patents FR-2,603,040 (US- 

5,463,894), FR-2,763,690, FR-2,772,477 filed by the appUcant are for example used. 

The sample initially saturated with a liquid A (brine for example) is placed in an 
etongate vessel or cup containing another fluid B of different density (oil for example). 
The vessel is festened to the end of a rotating arm and a centrifiigal force is applied 
20 thereto so as to study the displacements of the fluids in the sanq)le during at least two 
distinct stages. During a first drainage stage, the assenibly is subjected to a centrifi^al 
force exerted along the length of the vessel so as to exert thereon an e:q>ulsion force 
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BACKGROUND OF THE INVENTION 
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which tends to cause part of first fluid A to flow out. Simultaneously, part of fluid A 
penetrates the sample. The two fluids move in the sample to a position of equilibrium 
where the force due to the capillary pressure in the pores makes up for the centrifugal 
force exerted. A measuring sonde is placed in the vessel, on the side of the seunple. The 
sonde can be a capacitive type sonde for exanple, and it can comprise a metal rod 
insulated by means of a ceramic sheath. The capacitance between the rod and the 
conducting fluid (brine), which is proportional to the height immersed, is measured. 
With this system, the measuring accuracy is 1,5 % of the pore volume. The sonde 
detects the position of the interface between the two liquids in the vessel and transmits 
the measurement signals to a signal acquisition and control device comprising hydraulic 
liquid circulation means and an acquisition device. 

The various points of the capillary pressure curve (Pc) are obtained from 
measurement of the mean saturation of the sample at equilibrimn for given rotating 
speeds. Detemiination of the relative permeability curves (Kr) is achieved by means of 
the saturation variation with time between two (transient) speed stages. None of these 
curves is obtained directly ; the first one must be calculated and the second one requires 
a numerical flow simulator comprising (Fig.l 1) a simulation module which reproduces 
the effects of the experiment by means of a numerical model, and an inversion module 
allowing to determine capillary pressure and relative permeability curves by iterative 
adjustment of an initial niunerical model to the experimental results. 

An example of a weU-known flow simulator is the ACEX flow simulator based on 
Darcy's generalized equations, which is described in : 

- Chardaire-Riviere, C. ; 1992, In:^)roving the Centrifuge Technique by Measuring 
Local Saturations ; SPE 24882. 



1 

Although the centri&gation technique is widely used in the petroleum industry, 
because of its relative rapidity in relation to other existing methods, some 
implementation difficulties remain concerning determination of the capillary pressure 
curve as well as the relative permeability curves. 

5 For a given speed stage, the transient state is difficult to determine since a very slow 

production stabilization is observed. Search for this equilibrium leads to long-lasting 
operations spread over several weeks, whose implementation presents practical 
difficulties, considering especially the often high rotating speed of the centrifuge. 
Furthermore, the time required therefore is even longer if one tries to increase the pore 

10 volume of the san^le (by increasing the length thereof) in order to improve the 
saturation measurii^ accuracy. 

Inversion. i.e. adjustment of a numerical model to experimental results, is difficult 
to achieve because the information relating to the relative permeabilities is concentrated 
on the transient part of the production curve, Le. in a narrow range of esqwrimental 
1 5 values. The acceleration applied to the centrifuge has a preponderant effect on the form 
of this transient part, notably its duration. Parallel determination of the capillary 
pressure curve leads to applying several constant speed stages rather than a long 
acceleration, which reduces the transient curves. 

SUMMARY OF THE INVENTION 
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The modeUing method according to the invention aUows to obtain, with a 
considerable gain in time, the values of physical parameters of a geologic sample 
initially saturated with a fest fluid placed in a vessel and subjected to centrifugation in 
the presence of a second fluid. 
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It comprises measuring the amount of fluids displaced by centrifiigatioh and 
establishing satiiration curves, at fixed time intervals, and it is characterized in that it 
ccnq)rises, for each rotating speed stage, modelling of a saturation curve best adjusting, 
at each time interval, to the curve portion already obtained and corresponding to the 
various previous measurements, the value of the saturation at equilibrium being 
deduced from characteristics of the modelled curve. 

For each speed stage, an exponential type saturation curve is modelled by adjusting 
parameters, one of these parameters corresponding to the value of the saturation at 
equilibriimi. The saturation curve is modelled for example by means of a sum of 
e)q)onential fimctions depending on said adjustment parameters. 

Centrifugation speeds allowing to obtain an even distribution of points on each 
saturation curve are preferably selected. 

According to an embodiment, the method further conqirises : 

- using a flow simulator for modelling a relation between the relative permeability 
values in relation to the two fluids for a sample of determined dimensions and 
permeability, and the saturation values in relation to these fluids accordmg to various 
coeflHcients, and estabUshing a database coimecting parameters intended for adjustment 
to said coefficients, for several saturation values in relation to at least one of the fluids, 
and 

- determining relative permeabilities of the sample subjected to centrifiigation with 
reference to correspondence values extracted from the database. 




BRIEF DESCRIPTION OF THE DRAWINGS 

Other feature and advantages of the method according to the invention wUl be 
clear from reading the description hereafter of a non limitative example, with reference 
to the accompanying drawings wherein : 
5 - Figure 1 shows a centrifiigation device aUowing the method to be implemented, 

- Figure 2 shows an example of adjustment of a modeUing curve to an experimental 
production curve, 

. Figure 3 shows, by way of conqrarison, various modeUing curves obtained at 
different times of the same speed sts^e, 
10 - Figures 4a, 4b, 4c show the con?)ared evolutions respectively of a production curve, 
a variation curve of the time constant T2 involved in modelling of the simulation curve, 
and the variation curve of the final value Sa^ that can be predicted by simulation, 

- Figure 5 is a table showing an example of compared evolutions of several 
parameters durii^ the same ^>eed stage, 

15 - Figures 6a, 6b, 6c, 6d respectively show the evolutions in time of the saturation 
value of a sample, of coefficients Tl. T2 used as parameters in the simulated curve 
and the correspondii^ production curve, 

- Figure 7 is a table showing a second example of compared evolutions of several 
parameters during the same speed sts^e, 

20 - Figure 8a diagrammaticaUy shows a samration variation curve Sw and Figures 8b 
and 8c respectively show graphs counecti^g various coefficients and which correspond 
to two saturation variation zones. 



3NSDOCID: <GB_23S4690A ^l_> 



- Figure 9 shows the compared evolutions of saturation curves obtained according to 
the conventional method (Dl 1, D21) and according to the fast method according to the 
invention (D12, D22), 

- Figure 10 is a comparative table of numerical values for three speed stages, 
5 according to whether the conventional method or the accelerated method according to 

the invention is used, and 

- Figure 1 1 diagrammatically shows the organization of a flow simulator. 

DETAILED DESCRIPTION 
A centrifiigation device as described hereafter is used for implementing the method. 

10 It comprises (Fig.l) an elongate vessel or cup 13 intended for a sample. A 

capacitive type level-detectii^ sonde (not shown) is arranged in a tubular lateral 
chamber parallel to the direction of elongation of the sample and of the vessel. Any 
. variation in the level of the conducting fluid in the vessel leads to a capacitance 
variation of the sonde. Centrifuge II comprises a tank 9, an electric motor 10 whose 

15 shaft drives a hub 11 in rotation. Two identical arms 12 are mounted in opposition on 
hub 11. Two vessels or cups 13, 14 are mounted swivelling at the ends of arms 12 so as 
to spontaneously line up along the direction of the centrifugal force applied, and they 
balance each other in rotation. A porous rock sample to be measured, initially saturated 
with an electrically conducting fluid A, such as brine for example, is placed in vessel 

20 13. Vessel 13 is initially filled with an electrically insulating other fluid B such as oil. 

Two lines 15, 16 &stened along one of arms 12 connect vessel 13 to a rotary 
electro-hydraulic connector 17 carried by hub 11. First line 15 is used for injecting oil 
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into vessel 13. The second one, 16, is used for coUecting the brine drained out of the 
sample by the rotation. A cable 1 8 connects the measuring sonde to rotary connector 17. • 
By means of connector 17, the conductors of cable 18 are electrically connected to 
another cable 21. These two lines and this cable form a linking element 22 between 
5 vessel 13 and an external assembly E suited to coUect the signals coming from the 
sonde, to process the acquired data, to control fluid transfers and to actuate driving 
motor 10. 

The method conq>rises first selecting a certain number of rotating speed stages 
selected to obtain well-distra>uted points on the saturation curve. The operator can 
10 select them by referring to values obtained with previous tests for samples of the same 
type. Failing that, he can also quickly carry out a test on a sample by means of a Purcell 
type method well-known to speciaUsts, so as to establish a corresponding capillary 
pressure curve and, according to the form thereof to estimate the suitable speed stages. 

The fest access to the saturation vahies obtained with the method is based on 
15 analytic modelling of the fluid production stages of the sample subjected to 
centrifiigation. An exponential type modelling equation is used therefore. 

It has been found that most production curves can be modelled by the following 
equation : 

SiO ^sJl-iP- e-' ' ^' + (1 - P) • e-" )]+ 5^, 

where S., is the saturation at equilibrium, S^^ the initial saturation of the previous speed 
20 stage, t the time, P a weighting factor, and T„ Tj are time constants. 
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The lowest time constant Tl describes the transient part CI of the saturation curve 
(Fig.2), i.e. the dynamic effects. The highest time constant T2 describes stabilization 
zone C2, i.e. capillary effects. 

In order to estimate the saturation at equilibrium, long before the stage where its 
variations actually become lower than the measuring accuracy, coefficients Tl, T2 are 
sought which, applied in the equation above, generate a saturation curve that adjusts to 
the part of the production curve that is already known. 

At each new rotating speed stage, the control device performs a new measurement 
acquisition and a new pair of data is obtained (time, production). The programme then 
adjusts the bi-exponential mathematical modelling to the points curve already available 
(Fig.3). The evolution of time constant T2 is thus obtained, which determines the 
stabilization and final saturation Sa^, with time. In this programme, a maximum 
difference is fixed between two consecutive values of T2 and Se^, corresponding for 
example to the practical precision of the measuring device. When the calculated 
differences are below these values, it may be assumed that T2 and Ss,«, are sufficiently 
stabilized and therefore that the production curve is entirely determined, and the next 
speed stage can be dealt with. 

Their evolution is shown in Figs.4a to 4c. A very marked stabilization of time 
constant T2 and of the final saturation can be observed. This stabilization occurs as soon 
as the flow reaches the final plateau, i.e. as soon as the capillary forces become 
preponderant. Constant Tl stabilizes much fester since it only relates to the first part of 
the curve. 



The table in Figure 5 shows the evolution with time of four parameter sets P, Tl, 
T2, Sfiod calculated at four different times of the experiment. 

If the accuracy of the capacitive sonde measuring saturation S in the vessel 
containing the sample is of the order of 1.5 % of the pore volume for example, i.e. 
1 .5 % in saturation, it can be considered that the final saturation value is known after a 
3-hour experiment. This impUes a 2.6 % error on T2 and a 0.9 % error on S^, below 
the accuracy of the sonde. If one waited for effective stabilization of the saturation 
where the differences between successive measurements become smaller than the sonde 
accuracy, the speed stage should last about ten hours. The duration of the experiment is 
thus divided by three. 

The peak, then the stabilization observed on the simulated curves are found for 
parameters T2 and S«„^ from an experimental curve. The evolutions shown in Figs.6a to 
.6d are obtained. 

The table in Figure 7 shows similar results for another experimental case. 

After 3.5 hours, the final saturation value is known to 0.66 %, which is less than the 
measuring accuracy. If one had waited for stabilization, the experiment would have 
lasted for more than 7 hours. The duration of the test can thus be reduced by half, by a 
fector of the order of 2 or 3. It can thus be considered that, considering these two 
examples, when the value of T2 is stable to 5 % and the value of is stable to 1 %, 
the mathematical modelling parameters are sufficiently determined to allow the 
experiment to be stopped and the next speed stage to be dealt with. 
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Determination of the relative permeability curves 

Corey's equation, weU-known to specialists, is used for detennining the relative 
permeability curves for oil and water, Kr,nr and Kr,r respectively, from the respective 
saturations Snr and Sr obtained during the ejqjeriment : 

|^=5;-aHrf|^ = C' (2) 
where nl, .n2 are coefficients. 

A database is therefore formed, that can be presented in form of graphs for example, 
allowing to connect the coefficients of constant T2 to exponents nl, n2, for different 
values of these exponents. A sanq^le of determined length (6 cm for example) and with 
a given permeability (100 mD for example) is considered, and coefficients nl, n2 are 
varied between 1 and 5 each. 

Production curves S(t) are deterimned for various speed stages by means of the 
sunulator.As mentioned above, the values of P, Tl, T2 aUowing best adjustment of the 
modelled curve to the saturation curve variation are determined by applying relation 1. 
For all the value triplets (T2, nl, n2), value charts are formed or, more practically, 
graphs such as those shown in Figs.8a, 8b, correspondmg to different stages of each 
saturation curve (Fig.Sc). 

For a speed stage VI corresponding to a high water saturation Sw, it can be actuaUy 
observed that exponent n2 for oil has a predominant influence. The greater it is, i.e. the 
more the oil relative permeability is curved, the greater time constant T2. It can 
however be noted that, for this stage, the water exponent also plays a part, even though 
it is less marked. 



For speed stage V3 corresponding to a low water saturation Sw, the predominance 
of e^qjonent nl for water is very clear. E>q)onent n2 for oil has no influence any longer. 

For intermediate stages, it can be seen that the two exponents are by no means 
insignificant. At this stage, the flow is sensitive to the low relative permeability to each 
5 fluid. 

The tables or graphs formed can be used for samples of different leiigths and 
permeabilities, certain corrective factors being then applied- 
Using short samples allows to reduce the duration of a production stage. However, 
using a longer sample allows to spread the transient parts over the curves and thus to 
10 optimize determination of the relative permeabiUties. It is observed that, for samples of 
the same rock but of different lengths, a second degree polynomial type relation can be 
drawn between the stabilization time and the length. To go from a sample of length LI 
(the sample used for draNving the graphs) to a san^le of different length L2, the 
followii^ relation can be applied : 

15 A quasi-linear evolution of time constants Tl and T2 is observed as a function of 

the absolute permeabiUty of the sample. A corrective fector is therefore applied : 
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Validation of the method 



Two samples were centrifiiged, the jSrst one water wet, the second one rather oil 
wet, a first tune according to the conventional method, the second time according to the 
method described. 3 successive speed stages were applied : 470, 600 and 1050 rpm, 
5 determined after simulation tests, and the production curves shown in Figs.9a, 9b were 
obtained. 

Concerning equilibrium saturations, the results obtained are given in the table of 
Fig. 10. 

It can be noted that the final saturation accuracy greatly depends on the duration of 
10 the stage. 

The fest drainage for stages 1 and 3 thus leads to a final saturation that is slightly 
different fi-om the effective stabilization. However, the duration of these two stages is 
greatly reduced. Stage 3 on the other hand gives quite satisfactory results, but its 
duration is much less reduced. 

15 There is thus a compromise between measuring accuracy and diiration. After 

different tests, it is suggested to fix the variation criteria of parameters T2 and Sfi„ai to 
2 % and 0.5 % respectively. However, these values should be estimated according to 
each case. In fact, it can be noted that there are samples for which time constant T2 is 
difficult to stabilize when the final saturation is akeady well-determined. 

^0 Furthermore, the fact of not waiting for stabilization of the fluids before proceeding 

to a new stage does not disturb the shape of the production curve. In feet, a bi- 
exponential modelling can always be associated therewith. The precision of .the 
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modelling adjustment to the e3q)erimental curves remains very acceptable (of the order 
of 10-^). 

A modelling type involving exponential functions has been described. Other types 
of fimctions (splines for example) that can be connected to significant physical 
5 parameters of the san:5)le could however be used without departing from the scope of 
the invention. 
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CLAIMS 

1) A metiiod for feste- detenidim^ certain physical parameters of a geologic 
sanqjie initiaUy saturated with a first fluid placed in a vessel and subjected to 
centrifugation in the presence of a second fluid at speeds varying progressively in 
stages, con^ttising measuring the amount of fluids displaced by centrifugation and 
establishing saturation curves, characterized in that it conq>rises, for each rotating speed 
stage, modelling a saturation curve best conesponding at each tiine interval to the curve 
portion already obtained, correaponding to various measurements already performed, by 
selecting adjustment parameters, the value of each final saturation (Sen>i) being deduced 
from characteristics of the modelled curve. 

2) A method as claimed in claim 1, characterized hi that an ejqwnoitial type 
saturation curve is modeUed for each speed stage by adjustment of parameters (S.^, P, 
T), one of these parameters corresponding to the saturation value at equilibrium. 

3) A method as claimed in claun 2, characterized m that the saturation curve is 
modelled by means of a sum of exponential functioiis dependmg on said adjustment 
parameters (S,^ P, Tl, T2). 

4) A method as claimed in any one of the previous claims, characterized in that it 
comprises selectmg centrifugation speeds allowmg an even distribution of points to be 
obtained on each saturation curve. 

5) A method as claimed in any one of the previous claims, characterized in that it 
conpises: 
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- using a flow simulator for modeUing a relation between the values of relative 
permeability (K) respectively to the two fluids of a sample of determined dimensions 
and permeability, and the values of saturation (S) with these fluids according to various 
coefficients (nl, n2) and establishing a database connecting adjustment parameters (Tl. 
5 T2) to said coefficients, for several saturation values relative to at least one of the fluids, 
and 

. determining relative permeabilities (K) of Ibe sample subjected to centrifiigation 
with reference to correspondence values extracted from the database. 

6) A method substantially as hereinbefore described with reference to the 
1 0 accompanying drawings. 
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